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Abstract The analysis of 90° vortex-vortex scattering is extended to ^ 
scattering in all head-on collisions of n vortices in the Abelian Higgs model. 
A Cauchy problem with initial data that describe the scattering of n vortices 
is formulated. It is shown that this Cauchy problem has a unique global 
finite-energy solution. The symmetry of the solution and the form of the 
local analytic solution then show that - scattering is realised. 
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In recent years, analytic and numerical studies have uncovered quite un- 
expected scattering processes of soliton-like objects in (2 + 1) and (3 + 1) 
dimensions. The analytic studies were mainly based on the geodesic ap- 
proximation whose validity has recently been proved for vortices and 
monopoles ||, and were mainly concerned with the scattering of 2 objects. 
There are, however, some interesting analytic results for the scattering of 
more than 2 extended objects [|J || ||. Among the most interesting pro- 
cesses is 90° scattering of 2 vortices in a head-on collision. In Ref Jj]], existence 
theorems for this process were given, which act as an underpinning to nu- 
merical methods and approximation techniques. In this Letter this analysis 
is extended to - scattering in all head-on collisions of n vortices. 

The model we study is the Abelian Higgs model given by the Lagrangian 
density 

C = 1 -{D^){D^Y - \F, V F^ - \{\ $ | 2 -l) 2 . (1) 

<3> is the complex Higgs field, = d^Q — zA^Q, /i = 0, 1, 2, is the covariant 
derivative and the gauge fields are defined in terms of the real gauge 
potentials as F^ = d fl A u — d u A fl , fj,, v — 0, 1, 2. The indices are raised and 
lowered with the metric tensor g = diag(+l, — 1, — 1). The Euler-Lagrange 
equations are 



D^D^ +-<$>(] $ | 2 -1) = 0, <V^+ % -{<$>*{D v <$>) - = 0. (2) 



For all A the Euler-Lagrange equations have static, finite energy n-vortex 
solutions of the form ||, 

Ai(r, 9) = - eij x j n^-, A = 0, $(r, 6) = /(r) exp[zn0] (3) 



for i, j = 1, 2, where 

T 

dr \r dr 

,. } 

dr \ dr 



r^f-^V/ 2 (r)[a(r)-l] = 0, (4) 



2r^- (r^f) - 2n 2 f(r)[a(r) - l] 2 - Ar 2 f(r)[f 2 (r) - 1] = 0, (5) 



with 

/(0) = 0, a(0) = 0, lim f{r) = lim a(r) = 1. (6) 
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In the special case A = 1, it can be shown || that the solutions actually 
satisfy the first order Bogomolnyi equations, so / and a satisfy, 

r/'-n(l-o)/ = ^o / + / a -l = 0. (7) 

In this case, there also exists a 2n parameter family of static n-vortex so- 
lutions describing vortices located at arbitrary positions. The reason for its 
existence is the fact that for A = 1, the net force between static vortices is 
zero. 

We now study the scattering of vortices during the time from shortly 
before to shortly after the collision. The fields are taken to be of the form, 

$(t, x) = $(f) + t((x), A (t, x) = 0, Ai(t, x) = Ai(x) + tB t (x). (8) 

Here A) is taken to be the solution (H). It is assumed that equations 
(|2p can be linearized in (t(,tBi) for t G (— e, e), e < 1. These linearized 
equations are solved for 

C = nf(r)k(r) 

( Bl ,B 2 ) = (-n sin[(n-l)fl] ^ + ^ , -n cos[(n-l)fl] ^'^ + ^ ) 

(9) 

where k is the solution of the equation, 

r 2 k"(r) + rk'(r) - k(r) [n 2 + r 2 f 2 (r)] = 0, (10) 

with k(r)e r — > 1 as r — > oo. (The functions (||) have been used by Weinberg 
[ PH in his study of the zero modes of the static solutions.) 

We use the zeroes of |$| 2 to define the positions of the vortices. For the 
solution (|S|) one obtains 

| $ | 2 = / 2 (l + 2ntk cos(n9) + n 2 t 2 k 2 ) 

> f 2 (l-n | t | k) 2 

For t 7^ 0, | $ | 2 has exaclty n zeroes, namely at r = p and 9 = 0, — , . . . , 2 ^ n ~ 1 - >7r 
for t < 0, and at r = p and = -, ^, . . . , (2(ra ~ 1)+1K for t > 0. This 
shows that the solution (|9|) describes - scattering. Here p is the point where 
k(p) = An. That exactly one such p exists follows from the fact that k(r) is 
monotonic decreasing, and behaves like fc_ n r _?l for small r. 
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The problem with this approach is that the linearization, though plausible 
for times shortly before and shortly after the collision, has not been justified 
in an entirely rigorous fashion. Guided in our choice of initial data by the 
results just discussed we will now bring our discussion to mathematically 
rigorous conclusions. 

The first rigorous result concerns the global existence of the solution. We 
can now show that for certain initial data a unique global finite-energy time- 
dependent solution of equations (§) exists in the Lorentz gauge, d^A^ = 0. 
To do this, we first subtract a backround field ($, A^) and write, 



$(t, x) = + x), Ap(t, x) = A^x) + a^(t, x). 



(11) 



For the backround field we choose the static solution As initial data we 
choose: 



0(0, x) = 0, ao(0,x) = 0, aj(0,x)=0, fori = 1,2; 

d t (j)(0, x) = nf(r)k(r), d t a (0, x) = 0, 

[rk'(r) + nk(r)\ 



[rk'(r) + nk(r)} 



(12) 



d t ai(0,x) = — nsin[(n — 1)9}- 

d t a 2 (0,x) = — ncos[(n — 1)9] 

i 

For this choice of backround field and initial data, a unique global finite- 
energy solution exists, since all the conditions from Ref |]lTJ are satisfied. 

An essential element of the proof which is based on Segal's existence and 
uniqueness theorem fl2]l , is an iteration method. The method starts with 
writing the Cauchy problem in the form, 



d t ^ = -iAfy + J 

where ^* = {a ,p , ai,pi, a 2 ,p 2 , <j>, it*), = d a^ 
where the operator A is defined by, 



7T 



(13) 
and 



A 



( r o o o \ 

o r o o 

o o r o 

V o o o r j 



A 



m 



1 




(14) 
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The vector J can be calculated from equations (fj) ( and is given in pj ). 

The solution of the Cauchy problem (0)111) can be obtained as the solu- 
tion of the integro-differential equation, 

V(t,x) = e~ lAt ^{0,x) + C dsex.p\-iA(t- s)]J(^(s,x)). (15) 

Jo 

In turn this integro-differential equation is solved by the limit of a sequence 
of successive approximations \l/ n defined by the formula: 

m n+1 {t,x) = e~ lAt ^{0,x) + /* dsexp\-iA(t- s)]JC9Js,x)), (16) 

Jo 

where \l/o = ^(0,x), with the initial data flT2|). We now establish certain 
symmetries of the initial data \l/o an d use ([16|) to establish these symmetries 
for the successive approximations ^/ n , and finally for the solution of (|l3|). 

The first transformation we study is x — > x! = Sx where S is the orthog- 
onal matrix 

„ / cos^ -sin^ \ 

S = • 2™ 2n (17) 

I sin— cos— / v ' 

\ n n / 

Under this transformation the initial data change as follows: 

#(0,xO = M x *(0,f), 

with 

( I ? ) ,A = cos— /, 5 = sin— /. (18) 
V u 1 / n ft 

We can see that J(Mi9(s,x)) = MxJ(^(s,x)), [M 1} A] = and 

exp(-«At)Mi^ n (s,f) = Miexp(-zAt)* n (s, f), 

which implies that \l/ n (t,x') = M!\l/ n (t,x) for all n G A/". From this follows 
) = Mi^(t,f) for the solution ^. 
Next we study the reflection (£1,22) - ► (^1,-^2)- Under this transfor- 
mation the initial data change as follows: tp(t,xi, —X2) = M 2 ip(t, x%, Xq), 



Mi 



1 





V 



\ 

A -B 

B A 

1/ 



5 



where 



M 2 



/ -/ \ 
0-/00 

0/0 

V o o o c 1 



CV = V* 



(19) 



Furthermore, J (M 2 ^ (s , x)) = M 2 J(^f(s, x)). A gain we have [A/2, A] = 0, 
and ty n (t, Xi, — x 2 ) = M 2 ^/ n (t, xi, x 2 ). From this follows ^/(t, xi, — x 2 ) = 
M 2 ^ n (t,Xi,x 2 ), for the solution 

Under the transformations considered, all terms in the energy density: 



£ = - I A,$ | 2 +- I A* I 2 + z 4 + 2^ + g(l $ I 2 -l) 2 , (20) 



are invariant. This leads to the following conclusion: If by using functions like 
I $ \ 2 ,F^ or £ , there is a way of defining the positions (xf(t), x 2 (t)),a = 1,2, 
of exactly n separate vortices, these n positions must lie on n radial lines 
separated by an angle — with equal distance from the origin. (Below we 
will use the the minima of | $ | 2 to define these positions.) Furthermore, 
one of these radial lines must be the positive Xi-axis, or make an angle ^ 
with the positive Xi-axis. Any vortex that does not satisfy these conditions 
immediately leads to In — 1 other vortices, because of the symmetries of 
our solution. Since our solution is continuous, these positions will change 
continuously such that at t — the n positions coincide, and after the collision 
the vortices move again on the radial lines just described. Therefore, they 
can either go back on the radial lines they came in on, or go back on radial 
lines shifted by an angle -. We will study a further symmetry and use the 
Cauchy-Kowalewskyi theorem |HJ to show that the second case is realised. 

The last transformation we study is x — > Mx where M is the orthogonal 
matrix 



1 



1 



1 



1 



A 




(21) 



Under this transformation the initial data change as follows: 



#(0,M£) = M 3 $(0,f), 
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with 



M 3 






V 





c 

D 






-D 
C 










a 



C 
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cos —a , D 

n 



. vr 
sin —a. 

n 



-a J 



We can see that J(M 3 ^(s,x)) = -M 3 J(*(s, £)), {M 3 , A} 

exp(zit)M 3 = M 3 exp(-zAt). 



(22) 



and 



From this follows ty(—t,Mx) = M 3 fy(t,x), and we see that all terms in the 
energy density are invariant under the transformation (t,x) — > (— t, Mx). 
This establishes - scattering for n vortices. To show that, for small time t, 
|$| has exactly n minima (so that we can identify n vortices), we use the 
Cauchy-Kowalewskyi theorem . 

From Ref ||, we know that / starts with an r n term and a starts with 
an r 2 term. Equation (|T0| ) shows that asymptotically near r = 0, 

k ~ k_ n r- n + k n r n . (23) 

Using the same techniques as in Ref 0, we can now establish the analyt- 
icity of the initial data, and verify that all the conditions of the Cauchy- 



Kowalewskyi theorem |]T3[ are satisfied. We therefore have an analytic solu- 
tion near the origin. This solution leads to the following asymptotic expres- 
sion for I $ I 2 : 



[n/2] / 

|$| 2 (t,f) ~ f n {x\+xl) n +2ntflk. n E(-l)" +P 

p=0 \ 



n 
2p 



X 



n - 2p x 2 2 p +Kt 2 , 



(24) 



where we have assumed that x\ and xi are of order e, and t is of order e n , 
so that higher order terms can be ignored. ( K is a constant that can be 
expressed in terms of the coefficients of the expansions of a, /, and k. [n/2] is 
the highest integer that does not exceed n/2.) For t ^ 0, the approximation 
( |24|) of |$| 2 has exactly n minima which lie at r — (n \t\k_ n y/ n and 9 = 
0, 2*,..., for t < 0, and at r = (nltl^) 1 /™ and 9 = ^, (2(n ~ 1)+1)7r 

for t > 0. This establishes the ^ symmetry of the process. Our rigorous 
method is, however, not capable of following the minima for large time t. 
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